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In the present paper, we study a system of doped antiferromagnet in three dimensions at finite 
temperatures by using the t-J model, a canonical model of strongly-correlated electrons. We employ 
the slave-fermion representation of electrons in which an electron is described as a composite of 
a charged spinless holon and a chargeless spinon. We introduce two kinds of U(l) gauge fields 
on links as auxiliary fields, one describing resonating valence bonds of antiferromagnetic nearest- 
neighbor spin pairs and the other for nearest-neighbor hopping amplitudes of holons and spinous 
in the ferromagnetic channel. In order to perform numerical study of the system, we integrate 
out the fermionic holon field by using the hopping expansion in powers of the hopping amplitude, 
which is legitimate for the region in and near the insulating phase. The resultant effective model 
is described in terms of bosonic spinous and the two U(l) gauge fields, and a collective field for 
hole pairs. We study this model by means of Monte-Carlo simulations, calculating the specific 
heat, spin correlation functions, and instanton densities. We obtain a phase diagram in the hole 
concentration-temperature plane, which is in good agreement with that observed recently for clean 
and homogeneous underdoped samples. 
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I. INTRODUCTION 

Since the discovery of high-temperature superconduc- 
tors of cuprates, it has passed more than two decades l!]. 
Besides their high critical temperatures (T) of supercon- 
ducting (SC) phase transition, these cuprates have sev- 
eral interesting properties like anomalous properties in 
the metallic state, existence of Fermi arcs, etc[2|- To 
explain these properties, various theoretical approaches 
have been proposed Q. Although ample knowledge have 
been accumulated, we still do not have a theory that has 
been proved and accepted as a "right" one. 

The t-J model [J| is regarded as one of the canonical 
models for high-Tc cuprates. The model excludes doubly- 
occupied electron states at each site reflecting the strong 
correlations (Coulomb repulsion) among electrons. This 
condition makes it hard to get convincing and solid un- 
derstanding of the model such as its phase structure. 

The slave-particle approach[5| using the slave-fermion 
(SF) or the slave-boson representation has been proposed 
in order to treat this local constraint on the physical 
states faithfully. In the SF representation, each elec- 
tron is described as a composite of a charged spinless 
fermionic particle called holon and a neutral bosonic par- 
ticle with spin called spinon. The SF approach is known 
to be superior j6| (giving a lower ground-state energy in 
the mean-field theory) than its statistics-reversed assign- 
ment, the slave-boson representation, in the region with 
small hole concentrations S {6 is just the density of holons 
per site). 

We have studied the SF t-J model in path-integral 
formalism^. Let us summarize the results of Ref.f?!]. 



The local constraint is exactly respected by using the 
CP^ (complex projective) variables (we write it z^a in 
Sect. II) for spinous. The fermionic holons {4'x) are de- 
scribed by Grassmann numbers. In path-integral expres- 
sion of the partition function, fluctuations of variables 
along the imaginary time give rise to certain imaginary 
term in the action. By assuming the short-range anti- 
ferromagnetic (AF) order between the nearest-neighbor 
(NN) spin pair, we integrated over a half of the spinon 
variables, those sitting at the odd sites, assuming a short- 
range (SR) AF order to obtain an effective model. At 
the half filling {S — 0) the effective model reduces to the 
Heisenberg spin model. It favors the so-called resonating 
valence bonds (RVB), the NN spin singlet pairs with AF 
coupling. As holons are doped, the AF order are gradu- 
ally destroyed because hopping of holons is associated by 
hopping of spinous without spin flips, which breaks some 
RVB's. 

Also there arises an attractive force between NN holon 
pair reflecting the energy released by breaking RVB's. In 
fact, the NN holon pair breaks only 7(11) RVB's while 
a holon pair separated at longer distance breaks 8(12) 
RVB's in two(three)-dimensional lattice. In Ref.Q, we 
have introduced a hole-pair field, the condensation of 
which implies the SC state, and derived its Ginzburg- 
Landau (GL) model in the hopping expansion. At the 
mean-field level, this GL favors the so-called flux phase 
corresponding to (s -f id)-wave symmetry. 

The slave-particle approach intrinsically possesses 
U(l) gauge symmetry, since the electron operator is in- 
variant under the local and simultaneous rotation of 
phases of holon and spinon fields. The possible charge- 
spin separation phenomena[8| has a natural and poten- 



tially simple explanation such that the U(l) gauge dy- 
namics in cuprates is realized in the deconfinement phase. 
In fact, in the deconfinement phase, holons and spinous 
may appear as unbound quasiparticles moving indepen- 
dently due to the weak gauge force among them. 

The slave-particle approach has yet another advantage. 
The mean field theory based on the slave-particle repre- 
sentation is basically capable to describe various expected 
phases including the SC phase Q- However, the criticism 
to this result may be common to every mean field the- 
ory, i.e., the faithful evaluation of effects of fiuctuations 
around mean fields are missing. It is rather hard to eval- 
uate such effects analytically in nonperturbative manner 
because the model has local gauge symmetry as men- 
tioned and associated zero modes may give rise to strong 
effects in the infrared region. 

One may think that numerical studies may be one vi- 
able approach as the successful example of lattice gauge 
theory in high-energy physics demonstrates. However, 
straightforward numerical studies such as Monte-Carlo 
(MC) simulations of the SF t-J model in path-integral 
representation are still not feasible because of the notori- 
ous sign problem in the fermionic determinant generated 
upon integrating over holon variables. 

In the present paper, we shall revisit the t-J model on 
a three-dimensional (3D) lattice in the SF path-integral 
representation with the purpose to study its nonper- 
turbative aspects by numerical methods. In order to 
avoid the difficulty associated with fermionic determi- 
nant mentioned above, we derive an effective model by 
employing the hopping expansion to evaluate integrals 
over fermionic holons. It is an expansion in powers of 
the hopping amplitude of holons. An effective expan- 
sion parameter is i x (5, so the expansion is useful and 
legitimate at sufficiently low dopings 5. In the region of 
applicability of the hopping expansion, the wild behavior 
of fermionic determinant is suppressed in a natural way. 
We respect the structure of interaction terms generated 
by the hopping expansion, but consider their coefficients 
as independent parameters in a flexible manner. This is 
partly because these coefficients acquire renormalization 
via higher-order terms in the expansion. 

As mentioned, the action in path-integral representa- 
tion involves the imaginary part reflecting the imaginary- 
time dependence of the variables. This brings some com- 
plications in numerical approach. In Ref. [7] we have seen 
that the integration over odd-site spinon variables makes 
the spinon part of the resultant action real. In this pa- 
per, we avoid this imaginary part in another manner by 
simply considering the region of flnite T's; at sufficiently 
high T such that the dependence of variables on the imag- 
inary time may be neglected. We expect that each phase 
obtained at finite T's survives down to sufficiently low 
T's including T = 0, so the obtained phase diagram at 
finite T's is useful not only for itself but also for low T's 
down to T = 0. 

In the practical numerical study, knowledge and tech- 
niques developed in the study of lattice gauge theory of 



high-energy physics are helpful. By making MC simula- 
tions of the effective model, we obtain a phase diagram 
in the 5-T plane, which contains AF phase, SC phase, 
and metal-insulator (MI) transition. The overall phase 
structure is similar to that observed in experiments for 
lightly-doped materials |10|. 

The present paper is organized as follows. In Sect. II, 
we explain and set up the model in detail. The holon 
variables are analytically integrated out by means of the 
hopping expansion to obtain the effective model at small 
5''s and finite T's. The model includes several variables; 
(i) the spinon field Zxa^ (ii) the auxiliary field for spin- 
singlet (RVB) amplitude of NN spinon pair (we call it 
Ux^i in Sect. II), (iii) the auxiliary field for amplitude of 
holon and spinon hoppings in the ferromagnetic (FM) 
channel {Vx^), which works as an order parameter of the 
MI transition, and (v) the hole-pair field [Mx^) for su- 
perconductivity. We introduce the hole-pair fleld and 
include the associated GL terms to the effective action, 
respecting the NN attractive force between holons as dis- 
cussed in Ref. 7]. 

In Sect. Ill, we flrst study the case without the super- 
conducting channel (by neglecting the GL energy of hole- 
pair field) . We present the results of MC simulations for 
the corresponding model, which we call UV model. We 
calculated spin correlation function to study the AF tran- 
sition, and instanton densities of [/-field and y-field to 
study the decay of AF order and the MI transition. We 
locate the AF and MI phase transition lines. 

In Sect. IV, we study the full model including the SC 
channel and discuss SC phase transition together with 
AF and MI ones. We modify the coefficients of GL terms 
of hole pairs from the leading-order values of hopping 
expansion of Ref. Q so as to describe the d-wave SC ob- 
served in experiments instead of s + id one. This is be- 
cause the SC transition is expected (and actually verified 
later on) to occur in the metallic phase and the higher- 
order terms of the hopping expansion should be included. 
Our standpoint is that we regard the hole-pair part of the 
effective model in a flexible manner, i.e., its structure is 
suggested by hopping expansion but its coeflEicients are 
relaxed to study the region beyond the validity of the 
leading order of the hopping expansion. We find that 
the SC state occurs always in the metallic phase, whereas 
the AF long-range order (LRO) can coexist with the SC. 
There appear two phase transitions related with the SC. 
One is a primordial SC transition that stabilizes the am- 
plitude of hole pairs and gives rise to a pseudo-gap in 
holon excitation energy, while the other is a genuine SC 
transition reflecting a phase coherence of hole pairs asso- 
ciated with the Higgs mechanism. 

In Sect.V we present discussions and conclusions. We 
discuss that the present model offers us an interesting 
possibility of new description of a SC state in the frame- 
work of gauge theory with local interactions. 

In Appendix A we give some details of the hopping 
expansion of path-integral over holons. 



II. THE t-J MODEL IN THE SLAVE-FERMION 

REPRESENTATION AND HOLON HOPPING 

EXPANSION 

A. Path integral expression 

We start with the standard t-J model on a 3D cu- 
bic lattice[lli]. whose Hamiltonian is given in terms of 
electron operator Cxa at site x (= xi,X2,x^) and spin 
a [ = l(t),2(i)] as follows; 
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x^fi ' 'Jx . '^x'^x-\-fi 



(2.1) 



X,fl 



where 

^X(T ^ V ^xa^xa } ^xa-) 

*^^ ^ o X] ^l<7^<y<y'^^<7'^ (^ ■ Pauli matrices), 

a. a' 

fix = / ,C^^Cx(7- 



(2.2) 



/i(= 1,2,3) is the 3D direction index and also denotes 
the unit vector, a (1 = 2,2 = 1) denotes the oppo- 
site spin. The doubly occupied states (CJ^_|.CL |0)) are 
excluded from the physical states due to the strong on- 
site Coulomb repulsion. The operator Cxa respects this 
point. 

We adopt the slave-fermion representation of the elec- 
tron operator Cxa as a composite form, 



Cxa = -ipU 



(2.3) 



where ipx represents annihilation operator of the 
fermionic holon carrying the charge e and no spin and Uxa 
represents annihilation operator of the bosonic spinon 
carrying s = 1/2 spin and no charge. Physical states 
|Phys) satisfy the following constraint. 



[^al^axa + iplipxj Iphys) = |phys 



(2.4) 



In the salve-fermion representation, the Hamiltonian 
p.ip is given as 



J 

T 

X,fJ, 



y^ \{a)da)x+^j. ■ {a)Sa)x ~ {a^a)x+i_i{a)a)x , 
{a^a)x = '^a\^axa, {a' $a)x = '^a^xa^<ya'axa' ■ (2.5) 



We employ the path-integral expression for the partition 
function of the t-J model. 



Z = Trexp(-/3i?), /3 ; 



1 



feT 



(2.6) 



at finite T in the slave-fermion representation. This is 
done by introducing a complex number axaij) and a 
Grassmann number ipxi^) at each site x and the imagi- 
nary time T S [0,/3]. The constraint (|2.4p is solved f7!| by 
introducing CP^ spinon variable Zxai^), i.e., two complex 
numbers Zxi, Zx2 for each site x satisfying 



and writing 



(1 - tpxipx) ' Zxa 



= 1, 



1/2, 



(2.7) 



(2.8) 



It is easily verified that the constraint (|2.4p is satisfied 
by Eqs. (j2.8l ) and (12. 7p . Then, the partition function in 
the path-integral representation is given by an integral 
over the CP^ variables Zxa{T) and Grassmann numbers 

i'xij). 

We shall consider the system at finite and rela- 
tively high T's, such that the t— dependence of the 
variables Zxaji^x cire negligible (i.e., only their zero 
modes survive). Then the kinetic terms of Zxa-,'>px hke 
Zxdzx/dTjipxd'ipx/dT disappear, and the T-dependence 
may appear only as an overall factor /3, which may be ab- 
sorbed into the coefficients of the action and one may still 
deal with the 3D model instead of the four-dimensional 
model. Study of finite-T properties of the systems gives 
us an important insight into the low-T phase structure, 
for we can expect that ordered phase at finite-T generally 
survives at T = 0. 

In this way, the partition function Z of the 3D model 
at finite T's is given by the path integral^?,]. 



Z = I [dz] [di)] [dU] exp A, (2.9) 

[dz] = \{dzx, [dijj] = \{d,l:xdijx, [dU] =\{dUx 



'x^n 



X,fJ, 



with the following action A on the 3D lattice [13. [T3|. 

A^ Aaf + Ay^op + Asc, 



^AF " y XI {4+t,Ux,,zx + c.c. 



X,fl 



^hop — 'Z' / ^ \Zx-\-fiZxWxWx-\-fi + c.c. 



)-"l^/ 



a^,/i 



J/3 



^SC = -Y2^Px+t,Px\Zx+^,Zx\ 



(2.10) 



where 



Ux^ = expiiOxf,) e C/(l), 
Px = i>xipx, 

Zx+fiZx =^ Zx-{-f_i^iZxl "r Zx'^^^2Zx2 

^xl = Zx2, Z^2 = ~^£cl, 
Zx+fiZx = Zx-lrf^,2Zxl — Zx+fi,2Zxl- (2-11) 

The first term Ap^-p in the action A describes the AF 
coupling between NN spinous. We have introduced the 



U(l) gauge field Ux^ on the link {x,x + jj) as an auxil- 
iary field to make the action in a simpler form and the 
U(l) gauge invariance (explained below) manifest. The 
second term Ahop describes simultaneous NN hopping of 
a holon and a spinon keeping its spin orientation (i.e., in 
the FM channel). The third term Asc describes attrac- 
tive force between hole pairs, which we shall discuss in 
Sect.IID in detail. There are remaining terms 14|, which 
are irrelevant to discuss the global phase structure. 
The integration measure of z^a and Ux^i are 



//•OO /"OO 

dzx = 11 / dRezxa / dlvuzxa ■ S{^ 
„ J — oc -/ — OO 



^xa^: 



xa'^xc 



1), 



dlJx^, = 



^ de. 



Xfl 



27r 



(2.12) 



Grassmann variables f/'x anti-commute each other, 

['ipa;,^x']+ = [lpx,4'x']+ = [lpx,1px']+ =0. (2.13) 

The formulae of Grassmann integration [15J are 

diJxd^Px[l,i^x,^x,4'xiJx] = [0,0,0,1]. (2.14) 

The term m'^^ipxipx adjust the hole density to d as 
(t^x^x) = 6. (2.15) 

Therefore the parameter m works as (the minus of) the 
chemical potential. 

The action A is invariant under a local (x-dependent) 
U(l) gauge transformation with a gauge function AtJIGI]. 



^xa ^ 6 Zxai 




Ux^. ^ e-'^-+^Ux^,e 


-iX. 


ijjx -> e'^^^px- 





(2.16) 

B. AF and Ferromagnetic spinon amplitudes 

The gauge field Ux^ is related to the spinon field Zx as 
{Ux^.) ■ 



\^x+u^x\ 



(2.17) 



which is obtained by maximizing the action Axp. There- 
fore Uxfi describes the (c.c. of) phase factor of AF NN 
spin-pair amplitude z^j^^Zx of Ea. (|2.11[ ). In fact, one 
can integrate out Ux^ in Ea. (|2.9l) and obtain 



"dC/]exp(v4AF) = exp(icpi)) 



^cpi =^^oglo{ci\zl^^zx\), 

Xfl 



(2.18) 



where Iq is the modified Bessel function. The effective 
term A^pi should be compared with the original expres- 
sion Acpi of the GP^ model. 



Zcpi = / [dz]exp{Acpi), 
^cpi = -ir2^\^x+^,zx\ . 



(2.19) 



This GP^ model describes the t-J model without holes 
(c3 — 0, Asc = 0), i.e., AF Heisenberg spin model at 
finite T. Note that the amplitude z^^a^x between NN 
spinon pair reads explicitly as 



c+fj.,2Zxl — Zx+fj,^lZx2- 



(2.20) 



This expresses the amplitude of spin-singlet AF combi- 
nation of NN spinous, which is called the RVB. Both 
models with A^pi and A^pi have similar behavior and 
it is verified that they give rise to second-order transi- 
tions at certain ci and J[l7J. The parameters ci in the 
action (j2.10p are related to the original ones as|17| 



Cl 



J 13 for Cl > 1, 

(2J/3)i/2 for Cl < 1, 



(2.21) 



For the coupling C3, the relation is straightforward. 



C3 -- t/3. 



(2.22) 



Let us see the meaning of GP^ term Aaf (the AF spin 
coupling) and the hopping term ^hop (the i-term) fur- 
ther. For this purpose, it is convenient to introduce an 
0(3) spin vector field £x made of spinon Zx, 



L 



ZxCTZx 



/ ^ ^Xa^ (T(T' 
(7,0"' 



4-4 = 1. (2.23) 



The NN spin correlation 4+/^ -4 is expressed by the GP^ 
amplitudes (such as Zx+^Zx) as 



•4 = 21 



^X^fl^X I 



-1 



^ ^Fa:+/i^a; 



where we have used the identity, 



^x+^^x 



' \Zx+^iZx 



(2.24) 



(2.25) 



So if the spinon hopping amplitude Zx+^iZx which ap- 
pears in Ahop has an absolute value near its maxi- 
mum, \zx+ij.Zx\ ^ 1, then the NN spins are mostly FM 



'-x+fj. 






1. On the other hand, if the spinon RVB 



amplitude z* 



X-\-fI^ 



takes values with \z*^^Zx 



1, then 

the NN spins are mostly AF, 4-i-/^ -4 ^ —1- These two 
amplitudes satisfy the sum rule (|2.25p . The AF phase 
and FM phase are characterized by the LRO in the spin 
correlation function (4 '4)' ^'^^ they can coexist with 
each other as we see in the following sections. 



C. Holon hopping expansion and auxiliary field 14, 



X,fl. 



In Eq. (|2.9p . one can integrate out the fermionic holon 
field ipx by assuming small holon density i5 of Eq. (|2.15p . 
In this region, the hopping expansion of ^x is applicable 
as it is an expansion in powers of S. Some details of 



the integration of ipx are given in Appendix A. After 
integration over ip^ we obtain 



[d^p] exp(ylhop) = exp(^hop) 



7 ^ xf^'3\^Sr\- P 

^hop — I "o" / / y I ^x+fi^x I 



-H^ 



n( 



7 , _[_[ {Zx+fi--^x) 

x,fj.<i/ plaq. 



^2:+/j. ^xj I 



(2.26) 



The second term of Ahop denote the product of z^+^z^ on 
the hnk {x, x + fi) [z^z^+fi on the Hnk (a; + //, a;)] around 
the plaquette {x,x + iJ,,x + ^ + v^x + v) and the elhpsis 
denotes non-local higher-order terms. Both the first and 
the second terms favor FM couplings of NN spin pairs. 

Then we introduce a vector field Wx^ as an auxiliary 
field corresponding to ZxZx+^, 



{Wx 



\ZxZx-^fj,/ 1 



(2.27) 



by using Gaussian integration (Hubbard-Stratonovich 
transformation) as follows, 



exp(^hop) 



^--<y 



[dW]exp{Aw), 

2 



x,/j 



w. 



Xfll 



ZxWxu + C-C. 



' 7 ^ {zx+fi'^xyyxfi, 

X,^ 

*(?)' s: n 



w^ 



(2.28) 



x,fj,<i/ plaq. 



Estimation of the magnitude of Wxfj. is straightforward 
for T/J < 1 as 



[dz]e-i^\-^-^-^-^"\zx+^Zx\' 



1 



(2.29) 



Then we set 



Wx 



WVx 



w 



/jp' 



Vxfj. = exp{i(f>xf^) e U(l), 



(2.30) 



by ignoring the fluctuation of radial component of Wx 



and focusing on its phase, dWx 



dVx. 



So the correspondence (|2.27p becomes 



{V. 



x^l 



l^xZx-'rfi I 



d4>xf,/{2'K) 



(2.31) 



This simplification is based on the observation that the 
most relevant degrees of freedom in gauge theories are 
the phases of gauge fields on the links rather than the 
amplitude W because the latter has only massive exci- 
tations. The phase factor T4^ is a new U(l) gauge field 
that transforms under the gauge transformation (I2.16P 
as 



Vx^ 



JK 



y x^^ 



'i^x 



(2.32) 



Physical meaning of 14^ is obvious from the discus- 
sion given in Sect.IIB. It measures the phase part of 
SR FM spinon channel, the deviation from the AF or- 
der. Its coherent "condensation" induces coherent hop- 
ping of holons il)x in the FM spinon channel as the term 
Ahop shows, and therefore a MI transition into a metal- 
lic phase. More detailed discussion will be given in the 
following sections. 

The Ahop term is then rewritten effectively as follows. 



exp(JLhop) = / [dV]eyiY>{Av), 



^y = y ^ {Vxtj,Zx+^_,Zx + c.c.) 



X,fl 



+ y 2J {yxiyVx-\-„,fiVx+fi.,i,Vxfi-\-c.c.) 



XjfiKv 



M^]-n£^ 



(2.33) 



We have neglected the higher-order terms in Eq. (|2.26p as 
they have smaller coefficients for T/J < 1 with numeri- 
cal damping factors. However, effects of these non-local 
terms can be expected qualitatively. As they have all 
positive coefficients, all of them favor the order of the 
field Vxi_i and so the metallic phase. From this point of 
view, the critical hole concentration Sc of the MI tran- 
sition obtained by the numerical study in Sect. IV might 
give an overestimation for the true value. 

The parameters C4 and C5 in Ay are related to the 
original ones as 



C4 



C5 



J/3 



st^p 



6^ 



(2.34) 



In the following investigation of the phase diagram of the 
system, however, we treat C4 and C5 in more flexible man- 
ner as free parameters that are proportional to 6 and are 
increasing functions of /3 = 1/T. As most of phase tran- 
sitions in the present model appear in the region ci 3> 1, 
we identify T and S from Eqs. ([1211 12311) as 






(2.35) 

At this stage, the original partition function Z without 
Am is expressed as 

Z^Zuv^ f [dz][dU][dV] exp{Auv), 

Auv = Aaf{zx, Ux^,) + Av{zx, Vx^,). (2.36) 

This "UV" model describes the competition between 
the AF-RVB spin-pair amplitude Ux^ and the FM spin- 
hopping amplitude Vxfn the latter is generated by inte- 
gration over holon hopping. 



D. Hole-pair field Mxfi and the full model ^tuii 

As shown in the previous section, in the effective ac- 
tion (|2.10|) . there exists the term Asc that describes an 
attractive force between NN holes doped in AF magnets, 
or more precisely, in a SR AF background. This attrac- 
tive force comes from the J-terms in the Hamiltonian 
(|2.5p . Actually, the two holes with a mutual distance 
more than one lattice spacing break twelve AF bonds of 
spins, while a pair of holes at NN sites break just eleven 
AF bonds. Thus the NN hole pair is favored energeti- 
cally. To see it explicitly, we rewrite ^sc in Eq. ()2.10|) as 
follows, 



A 



sc 



•11 
2 



E 



■4'x 



-l-M 



^^x+fj.^x)Vx 



(2.37) 



Note that the holon-pair variable ipx+fii^x is accompa- 
nied with the RVB spinon-pair amplitude z* 



4-M 



This 



combination is nothing but Cx+fi,2Cxi — Cx+ij.,iCx2 in 
terms of electron operators. We expect that this attrac- 
tive force induces hole-pair condensation under certain 
conditions, and as a result, a SC state is generated. The 
main problem to be addressed here is whether the above 
attractive force is strong enough to generate a SC state 
in the region without AF LRO. 

In order to investigate a possible SC phase transition, 
we introduce a hole-pair field Mx^i as a complex auxiliary 
field describing the configuration of holon-pair accompa- 
nied with the RVB spinon pair at the sites x and x + fi. 
So Mxij, should satisfy 



(Mx^) - {^px+^{z*^^Zx)tpx)- 



(2.38) 



This hole-pair field M^^ is nothing but annihilation op- 
erator of spin-singlet electron pair sitting NN sites as 
mentioned. Explicitly, we use the Hubbard-Stratonovich 
transformation for Asc as Eq. (|2.28p . 



expl^ 



V': 



= / dMxf_cexp 



X + fJ-y^X+fi^xjV: 



MxuMx, 



+ 



— [Mx^il^x {zIj^^Zx) ^x+^, + ' 



(2.39) 



This assures us of Ea. (|2.38p . 

To study the effect of Asc, we start with Z ofEq. (P?TU) ) 
and rewrite Asc in the action by using Ea. (j2.39p . Then 
we integrate out the holon variables %px as in the previous 
case (without Asc there) to obtain the effective action, 
Afuii, where the suffix "full" implies Asc is taken into 
account. The partition function of the full model is now 
given as 

Zm, = Adz][d[/][dy][dM]exp(AMi), 

Afuii = Aery + Aa/ = Aaf + Ay + Am. (2.40) 



In addition to the action of the UV model of Eq. (12.361) , 
Afuii includes an extra term Am{zxtMx^i) that depends 
on Zx and Mx^i- 

We have calculated Am in the order up to 0((c3)'')[3]. 
Ea. (|2.39p shows that as ij^x and tpx+^. hop, they leave 
the factor Mx^z^j^ Zx-, i.e., Mx^i is always accompanied 
with the AF component of spinon, (c.c. of) z^j^^Zx- The 
hopping term Ahop itself supplies the FM component 
ZxZx+y{^ \zxZx+y\Vxy) aloug the link [x, x + v) ^x hops. 
In expressing Am we prefer to use Ux^i instead of z* 



X'\-fl'^ 



using Eq. (|2.17p . because it makes the gauge invariance 
of the system manifest. Then Mxfi appears in Am in 
the combination Mx^i,{z*^f^Zx) -- Mx^,\z*_^^^Zx\Uxf,- So 
we define a new variable. 



M,% = Mx^Ux^ 



i^x+^ii^x 



(2.41) 



and write Am in terms of M* and V:, 
Wc note that M. 



xfx , the latter is sup- 
plied by Ahop- Wc note that M*^ is not gauge- invariant 
and represents the "holon pair" at (x, x + fi) in contrast 
with gauge-invariant AIx^ for "hole pair" . 

In the practical calculations in Sect. IV, we focus on the 
phase degrees of freedom of M* , ignoring fiuctuations of 
the radial part of M*^ as in the case of Wxp. -> T4^ (the 
London limit). So wc set 

M*,.^Mexp{i^x^), 
M ~ \/holon-pair density -- S, (2.42) 

and dMxf_i = dipxf_i/{2Tr). We include M into the co- 
efficients of Am and treat Mr^ = exp(i(pxfi) as a U(l) 
variable. Furthermore, we regard \zxZx+ij.\ and l^^+^^a;! 
involved in Am as constants. They are also absorbed in 
the coefficients. The reason of this treatment is given 
below on the determination of the coefficients. 
Am is expressed as 



In terms of this Af * , 



Am = y E ^.'m^4*+.,^^. 

4 



x.^<v 



■*/Jl^ 



1,U ^X] 



VxvVx+v.,tiM*^^^^M*^ 



+Vxm:+,,^m:^^^,Vx^ 
+KuM:+,,^.Vx+f.,.Vx^ 

+M*Vx+u.^.Vx+^^uM* 



4 E M^^M^^^j'i: 



j_ M* 

-j-fl,!^ Xfl 



X^fl<l^ 



-l-c.c. 



(2.43) 



Each term in Am is schematically shown in Fig[T] The 
terms with the coefficients /i and /2 in Eq. (|2.43p describe 
the local hopping of the holon-pair field M*^, whereas the 
/3-term controls fluxes of M*^ penetrating each plaque- 
tte. These fluxes correspond to vortex excitations in the 
SC state. In other words, the /i and /2-terms induce a 
primordial SC state and a genuine SC state is generated 
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FIG. 1. Each term of Am of Eq.^^. The hnes with 
the reversed arrows indicate complex-conjugate variables, 
M^f^, Vx/j,- The gauge invariance under Ea. (|2.16p forces the 
arrows near each corner to make a divergenceless flow. 



by the /3-term. Numerical investigations in the following 
sections verify this qualitative expectation. 

As stated in Sect. I, we think that a SC state is to be 
realized in a metallic phase, i.e., beyond the region of ap- 
plicability of the leading order of the hopping expansion. 
So keeping the results of the hopping expansion for the 
three coefhcients fi of Aj^i is not suitable for discussing 
a SC state. For example, the coefficient /a is negative 
in the leading order of the hopping expansion, which fa- 
vors the s-l-id-wave SC. We examined higher-order terms 
of the hopping expansion and found that some of them 
generate a positive value of /a to support the d-wave SC 
as observed experimentally. So in the following numer- 
ical studies, we shall assume that /^'s are positive and 
proportional to S^, 



fuh,h^5^ 



(2.44) 



and treat their coefficients as positive free (phenomeno- 
logical) parameters 18]. In short, all the effect of |Afa;^t|, 
\zxZx+^L\ and \z*_^_^Zx\ in and near the SC state are in- 
cluded in the effective coefficients fi of Am- 

Also we have incorporated in Eg. (12.431 
structure of the 3D lattice of cuprates 



the layered- 
11| by introducing 



^/xz^ 



^Ufl 



{^ or v — 3) 



(2.45) 



The layered structure of the system is systematically in- 
corporated in the original Hamiltonian (12.11) by mak- 
ing the parameters t and J anisotropic. This induces 
anisotropics in the effective model that we have derived. 
Most of the terms are insensitive to the anisotropy except 
the /2-term, which is the reason that we treat Ahop and 
Ay in a symmetric manner. For the /2-term, the layered 
structure plays an important role to avoid frustrations 
and make the symmetry of SC to be d-x2_y2. 



III. PHASE STRUCTURE OF THE UV MODEL: 
AF AND MI TRANSITIONS 

In this section, we study the UV model with the action 
Auv = ^AF + Ay in Ea. (|2.36p by means of the MC 
simulations. The full model Af^w = Aaf + Av + Am shall 
be studied in Sect. IV. For MC simulations, we consider 
a 3D cubic lattice of the size V = L^ (L up to 30) with 
the periodic boundary condition. We used the standard 
Metropolis algorithm for local update. Average number 
of sweeps was 2 x 10^, and acceptance ratio was about 
~ 40%. 

To study the phase structure of the model, we mea- 
sured the internal energy E and the specific heat C, 
which are defined as 

E=-j^{A), C=j^{iA-{A)r), (3.1) 

as functions of the parameters ci, C4 and C5. We note 
that the C4 term and C5 term are related with each other 
because both are generated by C3 term. Below we respect 
this correlation by setting the parameter C5 as C5 a C4. 

By obtaining the locations of the phase transition lines 
by the peaks of C, etc., we get a phase diagram in the C4 — 
ci plane. Then we investigate spin correlation functions 
and instanton densities in order to identify the physical 
meaning and properties of each phase. To support this 
procedure, we also investigated fiuctuations of each term 
of the action by measuring the individual "specific heat" 
Cai defined by 



Ca, 



T.U. 



{A^)f 



1,4,5, 



Ai = A^-pi ^4.5 = C4 5— term in A\ 



(3.2) 



in Eq. (j2.43p the anisotropy parameter a^u, which is de- 



At C4 = (i.e., at a = 0), the system is reduced to 
the AF Heisenberg model with the action Aaf alone, 
which has a phase transition from the paramagnetic 
(PM) spin-disordered phase to the AF spin-ordered phase 
at ci ^ 2.8. We study how the location of this AF phase 
transition changes and whether new phases appear as the 
C4-term is turned on. It is naturally expected that the 
AF phase transition shifts to low-T region (large Ci re- 
gion) as the parameter C4 is increased because the C4 term 
favors FM NN spin coupling. 




FIG. 2. Total specific heat C and specific heat of each term 
CAi of Ea. (|3.2p as functions of C4 for ci = 3.5 and cs = C4/3.O. 
System size is L = 30. C has two peaks at C4 ~ 1.5, 3.0. Cai 
has a sharp peak at C4 ~ 1.5 suggesting the AF transition, and 
Ca4, Ca^ have peaks at C4 ~ 3.0 suggesting the MI transition. 



Let us first examine C and Ca^ as functions of ca for 
ci — 3.5. As we shall see later on, this value of ci belongs 
to relatively high-T region. In Figl2]we show the result 
for the case C5 = C4/3.O. We found no anomalous be- 
havior of E such as hysteresis, whereas C shown in Figj^j 
exhibits two sharp peaks at C4 ~ 1.5 and 3.0. We verified 
that each peak has a systematic system-size (L) depen- 
dence, so we concluded that both peaks show existence 
of second-order phase transitions. Cai of Figl2] exhibits 
a very sharp peak at C4 ~ 1.5. On the other hand, both 
Ca^ and Ca^^ exhibit a peak at C4 ~ 3.0. Then we con- 
clude that the AF phase transition takes place at C4 ~ 1.5 
and the MI transition at C4 ~ 3.0. 

The above conclusion may be confirmed by calculating 
the spin correlation function. In Fig|31 we show the cor- 
relation function Gs(r) of the 0(3) spin i^ of Eg. (12. 231) . 






(3.3) 



As we expected, at C4 = 0.7, Gs{r) exhibits an oscillatory 
behavior and has a staggered magnetization. 
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(3.4) 



So there exists an AF LRO at C4 — 0.7. This confirms 
that the phase transition at C4 ~ 1.5 is the AF transition. 
At C4 = 2.2 this AF order disappears and the system is 
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FIG. 3. Spin correlation function Gs{r) of Ea. (|3.3[) for ci = 
3.5, C5 = C4/3.O and L = 24. At a = 0.7, an AF LRO exists. 
At C4 = 2.2, the AF LRO disappears. At C4 = 3.8, a FM 
correlation appears as a result of existence of "free electrons" . 



in a magnetically disordered phase that we call param- 
agnetic (PM) phase. At C4 = 3.8, Gs{r) exhibits a LRO, 



Gs(r„iax)y^O, (FM phase), 



(3.5) 



which implies that the system is in the FM phase. So we 
obtain a picture of the phase structure for ci = 3.5 that, 
as C4 increases, the phase changes as AF -^ PM -^ FM. 
We also calculated instanton densities of the gauge 
fields Ux^ and Vx^- For example, C/- instanton density 
pu is defined for Uxfj — e'^=^" , O^fj. £ [— tt, tt] in the fol- 
lowing wavflTl [l9|. We first consider the magnetic flux 
©K/ii/ penetrating the plaquette {x,x + iJ,,x + i2 + v,x + iy), 
which is defined as 



^xaf 



(-4^ < e,^, < 47r). 



(3.6) 



Then we decompose Qx^l.v into its integer part Uxfiu, 
which represents the Dirac string (vortex line), and the 
remaining fractional part O^/ji/, 

^xiiu — '^T^nx^v -^Qx^I.v, (— TT < 82:^;^ < tt). (3.7) 

The [/-instanton density pu{x) at the cube around the 
site x-l-i-t-|-f| of the dual lattice is then defined as 



Pu 



(x) = ^^ X! '^A"'A(na;+/j,i/A - nx,,yx) 



^luX 



. / ^ ^iifXy^x+^.uX ^x^fX}-! 



(3.8) 



/i^i^.X 



where c^j/a is the totally antisymmetric tensor. From the 
above definition, we define the average instanton density 
Pu as 



Pu 



^E(i/'^(^)i)- 



L3 



(3.9) 



The I^- instanton density pv is defined similarly for Vx^- 



c 




11 

9 

7 

5 




4c4 5 



FIG. 4. Instanton densities pu and pv as functions of ca 
for Ci = 3.5, C5 = C4/3.O and L = 30. Arrows indicate the 
phase transition points determined by the specific heat (See 
Fig[2]). At the first phase transition point C4 ~ 1.5, pu starts 
to increase. On the other hand, at the second transition point 
Ci ~ 3.0, pv tends to vanish. 



The instanton density pu measures strength of fluc- 
tuations of the gauge field Uxp,- In the deconfinenient 
phase of Ux^i, fluctuations of Qx^v around its average 
©K/ii/ — are small and pu — Q. In the confinement phase 
of Ux^, on the other hand, Qx^lv fluctuates violently, and 
Pu has a finite value. Here we note that the confine- 
ment by Ux^ field gives rise to quasi-excitations that are 
gauge-invariant "composite particles" in the AF chan- 
nel. Such combinations include z*_^_ UxuZx, il'x+nUx^Zxa, 
etc. Similar interpretation holds for pv concerning to the 
gauge dynamics of Vx^. The confinement here works in 
the FM channel, and the possible gauge-invariant quasi- 
excitations are ^JxZxa = Cxcr, ^x^x, ZxaZxa' and their 
stretched versions such as ipx+fiVxfiZxcr, etc. 

In Fig|4] we show pu and pv for ci = 3.5. As we 
increases C4, pu starts to increase at the first phase tran- 
sition at Ci ~ 1.5. This result means the fluctuation 
of Uxii of AF NN spinon pairs become large, and the 
t/-conflnement spin-disordered phase appears. This re- 
sult is consistent with the interpretation based on Gg (r) 
above. On the other hand, at the second phase transi- 
tion at C4 ~ 3.0, Pv tends to vanish. So, for C4 < 3.0, the 
system stays in the y-conflnement phase and holons and 
anti-spinons are bound within electrons as ipxZxa- For 
C4 > 3.0, the system is in the y-deconfinement phase, 
and holons and spinous start to hop coherently and inde- 
pendently as low-energy excitations. This indicates that 
the phenomenon of charge-spin separation[8j takes place 
and also the system is metallic. 

Let us turn to the low-T region and see how the loca- 
tions of these AF and MI phase transitions change. In 
FigO we present the specific heat C and CAi for ci =6.5. 
We again found two peaks at C4 ~ 3.4 and 5.4. Figsl5)3, c 
show that both peaks develops systematically indicating 
that both phase transitions are of second order. Indi- 
vidual specific heat in FigjSJl shows that the peak of C 
at C4 ~ 3.4 corresponds to fluctuations of the C4 and C5- 
terms and so the MI transition, while the peak at C4 ~ 5.4 
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FIG. 5. Specific heat C as a function of C4 for c\ = 6.5 and 
C5 = C4/3.O. (a) C for L = 30. There are two peaks at 
C4 ~ 3.4, 5.4. (b,c) Each peak of C develops as the system 
size is increased. Results indicate that both phase transitions 
are of second order, (d) Specific heat Cai of each term. They 
indicate the transition at C4 ~ 3.4 is the MI one and the 
transition at C4 ~ 5.4 is the AF one. 



is generated by the ci-term and so the AF transition. 
Therefore the order of the AF and MI phase transitions 
along the C4 axis has been interchanged compared to the 
previous high-T case of C\ = 3.5. 

To verify the above observation, we calculated the spin 
correlations, Gs(r). The result is shown in FigjB] In the 
intermediate region 3.4 < C4 < 5.4, Gs{r) exhibits very 
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FIG. 6. Spin correlation function Gs{r) for ci = 6.5, cs — 
C4/3.O and L = 30. At C4 — 2.0, the oscillatory behavior 
around zero shows an AF order. At C4 = 4.0, there is an AF 
order in the FM background order. At C4 = 6.3, there is a 
FM order. 



interesting behavior, i.e., an AF correlation exists in a 
FM background. This implies coexistence of the FM and 
AF orders. We note that a coexistence of FM and AM 
orders was also observed previously in the CP^+Higgs 
boson model [201 ■ This model is a bosonic counterpart of 
the present model and the U(l) Higgs variable exp{iax) 
there plays the role of the fermionic holon variable ipx ■ 

In Fig|7] wc present instanton densities. Compared 
with the high-T result of FigHl the result again indi- 
cates that two phase transitions have interchanged their 
order along the C4 axis. As a result, there appears a 
range of C4 in which both pu and pv are small, which 
implies that spinous and holons hop here both in U and 
V channels. In other words, charges are transported by 
holons whereas spin degrees of freedom are transported 
by spinous both in the AF and FM channels. 
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FIG. 8. Phase diagram of the UV model for C5 — C4/3.O. (a) 
in the C4-C1 plane, and (b) in the 5-T plane. Each phase is 
separated by two transition lines; the AF transition line and 
MI transition line. All the transitions are of second order. 
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FIG. 7. Instanton densities pu and pv as functions of C4 
for ci = 6.5, C5 — C4/3.O and L = 30. Arrows indicate the 
phase transition points determined by the specific heat (See 
Fig[5]). At the first transition point C4 ~ 3.4, pv tends to 
vanish. On the other hand, at the second transition C4 ~ 5.4, 
Pu gets large values showing that the system enters into the 
[/-confinement phase. 



We repeated similar calculations for various values of 
ci and C4 and obtained the phase diagram of the model 
for C5 = C4/3.O. In FigslHl we present the phase diagram; 
FiglH^ in the C4 — ci plane and Fig|5}3 in the 6 — T plane, 
the latter is obtained from the former by using Eq. (12.351) . 

So far, we have studied the case of C5/C4 — 1/3.0. We 
also studied other values of the ratio C5/C4 and found 
similar phase diagram to that in Figs IS) As the value 
of C5/C4 is increased, the MI transition line shifts to the 
region of smaller S. This is expected from Ea. (|2.33p be- 
cause larger C5/C4 implies smaller critical value of C4, and 
therefore smaller critical S from Eq. (|2.35p . 

To estimate roughly the critical S of the MI transition, 
which we call <5mi, for real materials, one may put J ~ 
O.leV, i ~ 0.3eV so J/i ~ 1/3. Then we have 5 ~ (J/t)^- 
(C4/C1) ~ 0.1(c4/ci) from Eg. ip^Sl) . For ci ~ 10.0 (T ~ 
lOOK), the critical line of FiglS^ shows C4/C1 ^ 0.4, and 
this formula gives rise to Smi — 0.04. As discussed below 
Eq. (12.331) the higher-order terms in Eq. (|2.28p enhance the 
metallic phase, so the MI phase transition line is expected 
to be located in very underdoped region 5 <^ 1. 
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IV. 



PHASE STRUCTURE OF THE FULL 
MODEL: SC TRANSITION 



In this section, we study the fuh model of Eq. (|2.40p 
with the action Af^n = Aaf + ^v + ^m- Besides the 
AF and MI transitions observed in the previous section, 
we expect that the new term Ay^ in the action generates 
condensation of the hole-pair field M^,^ and/or the holon- 
pair field M*^ as the hole density (5 is increased. This 
condensation implies generation of a SC state. 

We studied the system Afuu by means of the MC simu- 
lations. As M*^ is a composite of holons at x and x + fi, 



we put /i,2,3 oc S' 



ex c| and M*^ 



€ U{1) as explained 



in Sect. II. Physically, the proportional constants fi/cf 
{i — 1,2,3) depend on the density of holes that actually 
participate in the SC fluid. We studied the system ^fun 
for various values of /i/c| and found that the system is 
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FIG. 10. Spin correlation functions for various values of C4 
for ci = 4.0, C5 = C4/3.O. L = 24. At C4 = 1.2 there is an 
AF order; At C4 = 2.4, no magnetic order; At C4 = 3.3, a tiny 
FM order; At C4 = 5.0, a FM order. 
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FIG. 9. Specific heat (a) C, (b) Ca^ , (c) Cf- as functions of 
C4 for ci — 4.0, C5 = C4/3.O and L = 24. Ca^h implies the 
fluctuatiuon of ^44-^5, and C/123 for A^f — Af-^ +Af^ +A/3. 



stable only for the case with small values of fi/cl. For 
example, the AF phase disappears at very small value of 
C4 for fi/cl ~ 0(1). In this section, we explicitly show 
the results for the case with fi — f2 — fz ^ 0.03 C4. 

Let us first study the high-T region first by choosing 
ci = 4.0, C5 = C4/3.O. In FigsEl we present various spe- 
cific heat as functions of C4. The total specific heat C in 
FiglS^ exhibits four peaks at C4 ~ 2.0, 3.2, 3.5 and 4.1. 
In order to identify the physical meaning of each peak, 
we show the individual specific heat C^^ in Figl9lD and 
C/; (i = 1, 2, 3) for the /^-term in Am defined similarly to 
Ea. (l3.2p in FiglS}:. From these results, it is expected that 
the first two peaks correspond to the AF transition at 
C4 ~ 2.0 and the MI transition at C4 ~ 3.2. The remain- 
ing two terms correspond to fluctuations of /-terms in 
the action, and therefore the SC phase transition. More 
precisely, the third peak at C4 ~ 3.5 in C corresponds to 
the /i, /2-terms and the fourth one at C4 ~ 4.1 to the 
/i, /3-terms. We shall comment on them later. 

In Figliniwe present the spin correlation function Gs{r) 
for various values of C4. It is obvious that only at C4 = 1.2 
the AF LRO exists. At C4 = 5.0, there is a solid FM 
order. This is consistent with the interpretation of four 
peaks above. 

In order to study the symmetry of SC state, we 
consider the quantity M2, the expectation value of 
M^fj^M^+^^ij {^^v), defined as 



M2 = -(M^+i,2M^i + Ma,+2.iM.^.+i,2 
o 

+ M^2M^+2,1 + Mx2M^i) + c.c. 



(4.1) 



In Fig[TT] we present M2. It takes negative values and 
starts to develop significantly at ci ~ 3.5, i.e., at the 
third peak of C. It is obvious that a d-wave correlation 
between adjacent hole-pair fields is generated beyond the 
third peak. 

We also measured the instanton densities pjj, pv and 
Pa/* • Pm* is defined in a similar manner to pu but by 
using the holon-pair field M*^{= M^f^Uxfj. -- 4'x+ti'4'x) 
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FIG. 11. Expectation value M2 of Eq. (|4.1|) for ci = 4.0, C5 = 
C4/3.O, L = 24. The result shows that d-wave correlation 
between adjacent hole-pair fields starts to appear at C4 ~ 3.5. 



instead of Ux^- Pm* reflects the vortex density of M* . 
These three instanton densities are shown in FigH^l By 
comparing FiglT2]with Fig|3]we see that the behavior of 
pu and pv is not influenced strongly by the existence 
of the /-terms, i.e., pu starts to increase at C4 ~ 2.0 
and Pv vanishes at C4 ~ 3.2. The M*-instanton density 
Pm* rapidly starts to decrease at C4 ~ 3.5 and vanishes 
at C4 ~ 4.1. We think that the SC phase transition, 
which is signaled by vanishingly small pM* , takes place 
at C4 ~ 4.1. 

From the above numerical calculations, we understand 
the physical meanings of the two peaks at C4 ~ 3.5, 4.1 
as follows. The third peak at C4 ~ 3.5 are generated 
by the /i and /2-terms and is located just after the MI 
transition at C4 ~ 3.2. After the MI transition, the holon- 
hopping amplitude Vx^ becomes stable, and as a result, 
these /i and /2 terms start to correlate the phases of a 
pair of adjacent link fields M* . In fact, these two terms 
in Eq. (|2.43p need a stabilized Vxf^ to let M* stabilize. 
FigfT2l shows that pM* at C4 ~ 3.5 is still large. So this 
effect is not strong enough to stabilize the holon-pair field 




M*^ completely at this region of C4. In order to suppress 
vortex excitations of the holon-pair field M* (making 
Pm* small enough), sufficient amount of the /3-term is 
necessary. The fourth peak of C at C4 ~ 4.1 corresponds 
to the critical value of /s to realize such M* stabiliza- 
tion with phase coherence and generation of SC. These 
consideration leads to our conclusion that the genuine SC 
starts at the fourth peak C4 ~ 4.1. 

From the above consideration, we expect that the re- 
gion between the third and fourth peaks corresponds to 
a primordial SC state. In this region we expect that 
holons acquire a pseudo-gap. In fact, as M*^ couples to 
Tpx as M*^'ipx+tJ.ipx, finite expectation value of M*^ sup- 
plies fermion-number nonconserving hopping processes 
effectively. Together with the fermion-number preserv- 
ing hopping term supplied in Ahop these processes give 
rise to a gap in excitation energy of holons [2l|. 

Furthermore, from the local gauge symmetry of the 
system, the terms like M* Zx+fj.Zx are also to be gener- 
ated by the renormalization effect of high-energy modes 
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FIG. 12. Instanton densities pu, pv and pu* as a function 
of C4 for c\ = 4.0, C5 = C4/3.O and L — 12. Arrows indicate 
the locations of four peaks in C of Fig[9] 



FIG. 13. Specific heat C, Ca^ and C/^ as functions of C4 for 
c\ — 6.5, C5 — C4/3.O and L — 24. 
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FIG. 14. Spin correlation functions G3{r) for various values 
of C4 with ci — 6.5, C5 = C4/3.O. At 04 = 2.7 there is an 
AF order; at C4 = 4.5, an AF order in a FM background; at 
C4 = 6.5, and C4 — 8.0, a FM order. 



of Zx and ipx- Then the spinon field z^ also acquires 
an extra contribution to its pseudo-gap, irrespective of a 
possible pseudo-gap expected by the mixing of two chan- 
nels, ciUxf_tZ*^^Zx and C4T4^Za:+^Za:- Anyway, the phys- 
ical properties of that state such as excitation spectrum 
is interesting and should be reserved as a future problem. 
Next, let us study the system Afun at lower-T region by 
setting ci = 6.5, C5 = C4/3.O. Behavior of various specific 
heats, C, Ca^, C/. are shown in FiglT31 There are again 
four peaks in the total specific heat C at C4 ~ 3.4, 3.6, 4.0 
and 5.6. From the behavior of Ca^ and C/. the first peak 
at C4 ~ 3.4 corresponds to the MI transition, the peak(s) 
at C4 ~ 4.0(and 3.6) to the SC transition, and the fourth 
peak at C4 ~ 5.6 to the AF transition. The order of these 
transitions is different from that at the previous high- 
T case as we have already seen in the UV model. In 
order to verify the above identification, we calculated the 
spin correlation functions Gs{r), the expectation value of 
adjacent hole-pair field M2, and the instanton densities 
as before. We show the results in FigsfHl [T51 and [T51 





4 6 C4 



FIG. 16. Instanton densities pu, pv and pm* as functions of 
C4 for ci — 6.5, C5 = C4/3.O and L — 12. Arrows indicate the 
phase transition points determined by the specific heat (See 
Fig[9|. Their behavior is consistent with the phase transition 
discussed in the text based on Fig[9l 



These results support the interpretation of each phase 
given above. 

In FigfTTl we present the obtained phase diagram of 
the fuU model Zf^ii of Eq. (P^ in the 6-T plane. Each 
phase are separated by three transition lines for AF, MI, 
and SC transitions. The SC phase always exists inside 
the metallic phase, whereas there is the coexisting phase 
of the AF and SC at the low-T region. In addition to 
these three lines, one may add the line corresponding to 
the primordial SC transition as the line of pseudo-gap 
generation. Except for the pseudo-gap transition, which 
seems not to be a sharp transition in experiments, this 
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FIG. 15. Expectation value AI2 for ci = 6.5, C5 = C4/3.O, and 
L = 24. The result shows that d-wave correlation between 
adjacent hole-pair fields appears at C4 ~ 3.6. It is interesting 
to observe that the correlation decreases slightly in the region 
without the AF LRO C4 > 5.6. 



FIG. 17. Phase diagram of the full model of Eq. !\2A0\i in the 
S-T plane. The three lines for AF transition, MI transition 
and SC transition separate each phase. All the transitions are 
of second order. One may add the pseudo-gap transition line 
as the fourth line. 
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phase diagram is consistent with that observed experi- 
mentally for homogeneous clean underdoped samples [lO|. 



V. CONCLUSION AND DISCUSSION 

In the present paper, we have studied the phase struc- 
ture in the underdoped region of the t-J model by using 
the slave-fermion representation. In this formalism, the 
AF-insulator phase naturally appears and it is expected 
that beyond a critical hole concentration 5mi the coher- 
ent hopping of holes is generated and the system enters 
into the metallic phase. This phenomenon was previ- 
ously studied by the mean-field theory, and the critical 
hole concentration was estimated as 5mi — [23 . 

We investigated the system by integrating out the 
fermionic holon field by the hopping expansion, which 
is legitimated for the region in and near the insulating 
phase, and then numerically studied the AF, MI and SC 
phase transitions. The obtained phase diagram is con- 
sistent with that observed experimentally for clean and 
homogeneous samples at small hole concentrations. The 
present study also implies that the observed pseudo-gap 
corresponds to a primordial formation of SC order pa- 
rameter Mxi_i. 

For the SC phase transition, we have treated the coef- 
ficients of effective action in more flexible manner than 
the original hopping expansion although we maintain the 
structure of interaction terms. As explained, this is be- 
cause these coefhcients certainly acquire renormalization 
and even change their signature as we go into the SC 
state. Some of our results in the present paper may 
reflect this flexibility, i.e., they may not be possible in 
the original t-J model due to the restrictions among the 
coefficients. The pseudo-gap transition might disappear 
(merge to the genuine SC transition) with different treat- 
ments of the coefficients. Even in such case, the results 
obtained in the present paper have important meaning 
as the knowledge of a reference system to the t-J model 
and other canonical models of the high- Tc materials. 

Concerning to the SC order parameter, we proposed 
gauge-invariant M^^, for hole pairs as a most direct 
possibility i7l]. We have calculated its correlation function 
(MxfiMy^), but found no LRO of M^fj, even in the SC 
phase. We understand this in the following way. If one 
could calculate this correlation of the t-J model exactly, 
one would have LRO in the SC state. The effective model 
in exact treatment certainly contains a lot of nonlocal in- 
teraction terms among M^^,, although their coefficients 
are small. Nonvanishing LRO is to be supported by these 
nonlocal interactions. However, the present model trun- 
cates the effective interaction terms to short-range ones, 
and so fails to produce LRO of M^fj- 

However, the study of lattice gauge theory [23j provides 
us with a viable alternative of describing a SC state. 
An effective system may involve only short-range inter- 
actions but it may generates the Higgs phase in which 
Meissner effects takes place actually. The price to pay is 



that there are no local order parameters to signal LRO. 
Our present model with the action Am is just a such 
model. Because our gauge-noninvariant M*^ for holon 
pairs has vanishing correlations due to gauge-invariant 
action due to Elitzur's theorem[2J], one need to introduce 
complicated nonlocal order parameters [2^ to show that 
some kind of LRO exists. Here we note that existence 
of LRO is a beautiful theoretical criteria to demonstrate 
SC phenomenon, but not a necessary condition. A simple 
and direct proof of a SC state may be to measure the mass 
of external electromagnetic field and demonstrate the 
Meissner effect, i.e., the Higgs mechanism. We have not 
made such a proof, but the existence of anomalous peak 
of the specific heat certainly demonstrates a new phase, 
which should corresponds to the Higgs phase. In fact, 
we have considered a U(l) Ginzburg-Landau modelt2(|, 
which is obtained from Am of Eq. (I2.43[) by putting T4^ to 
a certain constant. So the model loses gauge symmetry 
or viewed as a gauge-fixed version. The MC simulation of 
this model certainly exhibits a Higgs phase for sufficiently 
large fi in which the correlation functions {M*^M*^) ex- 
hibit a LRO. Let us summarize the situation. Because the 
faithful effective model of M^fj. is full of nonlocal interac- 
tions, we replace it by a short-range model. By sacrificing 
the LRO of gauge-invariant local order parameter, we are 
able to obtain the new phase. The analysis of the related 
model and the experience of lattice gauge theory strongly 
indicate this phase is a Higgs phase which is necessary to 
support SC. 

The reason why we integrate out the fermionic holon 
field analytically is obvious, i.e. it is technically diffi- 
cult to study fermion systems by numerical methods. In 
recent years, however, it has become possible to numer- 
ically simulate relativistic fermion systems and therefore 
it is important and also interesting to study the MI phase 
transition in the present system by means of those sim- 
ulation methods. This problem is under study and we 
hope that the result will be reported in a future publica- 
tion. Even in such a situation, the content of the present 
paper may be useful as some basis and a reference to 
obtain further understanding of physics of high- Tc su- 
perconductors. 



Acknowledgment 

This work was partially supported by Grant-in- Aid for 
Scientific Research from Japan Society for the Promotion 
of Science under Grant No.20540264. 



Appendix A: Holon-fleld integration 

In this appendix, we show some details of holon-fleld 
integration to derive Eg. (12. 261) . The same techniques 
are apphcable to derive Am in Eq. (|2.43p . It is use- 
ful to start with the original path-integral expression|7| 
in which Grassmann number ^xiT) is a function of the 



imaginary-time t. This is because the ordering of vari- 
ables is crucial to obtain the correct results. Then the 
relevant integration reads as 
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where we have used the following Green function of the 
hopping expansion, 



{lpx{Tl)iJx{T2)) 
g-m(ri-r2) 
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V2-ri)]. (A.2) 



In Eq. (jA.2l) . m is the chemical potential and there holds 
the relation, 



5={Mr + 0)Mr)) 
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(A.3) 
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